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Three  aspects  of  wave  propagation  in  a  random  medium  are  considered.  The  first  is  the 
scattering  of  a  plane  compressional  wave  by  a  spherical  inhomogeneity  of  specified  radius.  This 
result  is  then  averaged  over  a  Rayleigh  distribution  of  sphere  radii.  As  a  second  topic,  the 
scattering  by  a  medium  in  which  an  elastic  parameter  varies  randomly  is  considered.  When  the 
randomness  is  assumed  to  be  described  by  a  Gaussian  con-elation  function,  it  is  found  that  the 
frequency  and  angular  dependence  of  the  scattering  is  the  same  as  that  obtained  in  the  previous 
calculation  with  the  Rayleigh  distribution  when  the  average  diameter  of  the  spheres  is  set  equal  to; 
the  con-elation  distance  appealing  in  the  random  medium. 

The  third  topic  considered  is  the  effect  of  the  randomness  on  the  focusing  of  the  wave.  The 
total  incident  sound  field  may  be  viewed  as  being  composed  of  a  sum  of  a  coherent  field  reflected 
from  an  object  of  interest  and  an  incoherent  field  resulting  from  the  scattering  due  to 
inhomogeneities  in  the  surrounding  medium.  The  fluctuations  in  the  imaging  of  the  incoherent 
field  can  be  decreased  by  increasing  the  size  of  the  receiving  system  and  this  dependence  on 
receiver  size  is  determined  for  a  square  receiver.  It  is  found  that  most  of  the  advantage  to 
increasing  receiver  size  will  have  been  achieved  when  the  receiver  has  been  increased  to  a  size 
equal  to  about  twice  the  correlation  length  of  the  fluctuations  in  the  incident  field. 


A.  SCATTERING  BY  A  SPHERICAL  INHOMOGENEITY  WITH  AVERAGING  OVER  A  RAYLEIGH 
DISTRIBUTION  OF  SCATTERING  RADII 

The  sound  pressure  field  scattered  by  a  spherical  inhomogeneity  in  a  solid  medium  is 
examined.  The  wave  motion  is  described  by  a  standard  theoretical  model  used  in  geophysical 
prospecting'1)  in  which  shear  waves  are  ignored.  As  is  well  known,  the  angular  distribution  of  a 
scattered  field  is  sharply  peaked  when  the  circumference  of  the  scatterer  is  greater  than  the 
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wavelength  of  the  incident  wave.  When  the  average  scattering  produced  by  an  assemblage  of 
scatterers  of  different  radii  is  considered,  however,  the  variation  in  the  average  field  can  be 
expected  to  be  smoothed  out.  For  wavelengths  larger  than  or  comparable  to  the  circumference, 
there  is  little  angular  dependence  and  hence  little  effect  due  to  averaging.  The  averaging  is 
performed  here  with  respect  to  a  Rayleigh  distribution.  In  the  second  part  of  this  report  it  is  shown 
that  the  average  intensity  scattered  by  a  random  medium  with  a  Gaussian  correlation  function  has 
the  same  dependence  on  frequency  and  scattering  angle  as  that  obtained  for  the  spherical  scatterers 
when  the  correlation  length  is  set  equal  to  the  rms  diameter  of  the  spheres. 

A.l  The  Theoretical  Model 

As  noted  above,  we  use  a  simple  model  in  which  shear  waves  are  ignored  and  only 
compressional  waves  are  considered.  The  governing  equations  are  taken  to  be 


dt2 

P  =  -BV 


p--  =  -VF 


(1) 


where  P(r,t)  and  £(r,t)  are  the  sound  pressure  and  elastic  displacement  in  the  medium, 
respectively.  The  parameter  p  is  the  density  of  the  medium  and  in  order  to  keep  the  analytical 
model  as  simple  as  possible,  is  assumed  to  be  a  constant.  The  bulk  modulus  B  is  written  in  the 
form  B  =  Bo  -  8B  where  Bo  is  the  constant  background  modulus  and  the  scattering 
inhomogeneity  is  taken  to  be 


8B  = 


r  <  a 
r>  a 


(2) 


where  8Bo  is  a  constant. 

The  wave  equation  that  now  follows  from  Eq.(l)  has  the  form(2) 


V2P 


pd2P 
dt 2 


(3) 


Confining  attention  to  a  wave  at  a  single  frequency  (0  by  writing  P(r,t)  =  p(r)e-‘wt ,  we  find  that 
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Eq.(3)  has  the  form 


(4) 


where  Co2  =  Bo/po  • 

The  scattered  field  due  to  an  incident  plane  wave  pjnc.  =  poeik*r  is  given  by 


(5) 


where  the  integradon  is  performed  over  the  spherical  volume  of  the  scatterer.  At  distances  that  are 
a  number  of  wavelengths  away  from  the  scattering  volume,  the  expression  for  the  scattered  field 
reduces  to 


Psc 


— f  d3r'e~iks'r  p(r') 
4k  B()  r  iv 


where  ks  is  in  the  direction  in  which  the  scattering  is  observed  and  has  a  magnitude  equal  to  k  = 
cd/co  .  For  weak  inhomogenieties,  one  may  expect  the  field  within  the  scatterer  to  be  nearly  equal 
to  the  incident  field  and  thus  set  p(r')  =  pinc.  =  poeik*r'  in  the  integrand  of  Eq.(6).  The  scattered 
field  is  then 


Psc  _  ^  SB,  e  r  (pr'e~' K*r' 
Po  4tc  B0  r  Jv 


(7) 


where  K  =  ks  -  k  .  The  integral  is  now  readily  evaluated  and  yields 


P*L  =  k 2®q.£^  ^3  h(Ka) 

Po  Bfl  r  I  Ka  _ 


(8) 


where  ji(x)  is  a  spherical  Bessel  function  given  by  ji(x)  =  sin  x/x2  -  cos  x/x  and  K  = 
2ksin  (0/2)  in  which  0  is  the  scattering  angle,  i.e.,  the  angle  between  ks  and  k  .  Polar  plots  of 
this  result  for  various  values  of  ka  are  shown  in  Fig.  1.  As  expected,  side  lobes  in  the  scattered 
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field  do  not  appear  until  ka  =  2k /A,  becomes  greater  than  unity.  Since  interest  in  the  present 
project  is  confined  to  low  frequencies,  plots  of  higher  ka  values  have  not  been  considered. 


A.2  Averaging  over  a  Rayleigh  distribution  of  scattering  radii 

As  seen  from  Fig.  1,  the  angular  dependence  of  the  scattered  wave  becomes  quite 
pronounced  when  ka  >  1  .  The  net  effect  of  the  scattering  due  to  an  assemblage  of  scatterers  of 
different  sizes  should  be  a  smearing  out  of  these  sharp  peaks.  When  the  radii  are  assumed  to  vary 
according  to  a  Rayleigh  distribution,  defined  by 


r  /  \  _  ^  —a2  12 

/r  (^0  ~  2^ 

dr\ 


a0 


(9) 


where  ao  is  the  rms  size  of  the  scatterers,  we  obtain  the  average  of  the  scattered  field  given  in 
Eq.(8)  in  the  form 


■Bo  r 


(10) 

/ 


where 

(ll) 

When  ji  is  written  in  terms  of  an  ordinary  cylindrical  Bessel  function  as  ji(x)  =  (7i/2x)1/2J3/2(x), 
the  integral  in  Eq.(ll)  takes  on  a  standard  tabulated  form!3)  and  yields  the  much  smoother 
distribudon 


/ik\  =  (12) 

\PoLayl  V2  B0  hr  ' 

which  is  shown  in  Fig.  2.  Since  0  <  0  <  n  ,  no  sidelobes  are  generated  by  the  term  sin  (9/2)  in 
this  result.  For  ka  <  1  ,  the  averaging  is  seen  to  produce  little  effect  while  for  larger  ka  values 
the  peaks  are  smoothed  out. 
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Note  that  the  squared  magnitude  of  the  scattered  field  is  given  by 


(~) 

2 

n 

\  P®  /  Rayl. 

~2 

\  B0  V 

Of 


(13) 


In  the  next  section  it  will  be  shown  that  this  same  dependence  on  wave  number  k  and  scattering 
angle  0  is  obtained  for  the  scattering  by  a  random  medium  when  the  randomness  is  described  by  a 
Gaussian  correlation  function. 


B.  SCATTERING  BY  RANDOM  FLUCTUATIONS  IN  THE  BULK  MODULUS 


We  now  consider  the  scattering  per  unit  volume  due  to  random  fluctuations  in  the  bulk 
modulus  of  the  medium.  Returning  to  Eq.  (4)  and  replacing  5Bo/Bo  by  a  randomly  varying 
function  f(r),  Eq.  (5)  is  replaced  by 


i*|r-r'| 


r-) 


(1) 


where  the  integration  is  now  over  the  entire  random  medium.  Again  using  the  far  field  assumption 
r  »  r1,  we  have 


/Lc 

Po 


,  2  ikr 

k  e 


d:'r'f(r')e~iK’r 
4k  r  iv  V  ’ 


(2) 


The  scattered  energy  is  proportional  to  Psc(r)*Psc(r)  =lpsc(r)i2 .  We  now  treat  the  scattered  wave 
as  a  random  variable  and  form  a  statistical  average  that  is  proportional  to  the  scattered  energy  by 
considering 


Po 


(3) 


Setting  <f*(r')f(r")>  =  |i(|r'  -  r"l),  i.  e.,  assuming  that  the  correlation  of  fluctuations  in  the 
random  variation  is  isotropic  and  homogeneous,  we  can  introduce  the  coordinates  R  =  (r'  +  r")/2 
and  p  =  (r*  -  r")  to  obtain 
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-1vJ/VkM|p» 


(4) 


where  V  is  the  total  volume  of  the  scattering  region.  Carrying  out  the  integration  over  angles,  we 
obtain  the  scattering  per  unit  volume  in  the  form 


(Psc 

’}_  I 

f— ] 

2 

47t  r 

Vpl  1 

k4kj 

rK* 

pdpp(p)sin  Kp 


(5) 


The  remaining  integral  may  now  be  evaluated  for  various  assumed  forms  for  the  correlation 
function  |i(p)  .  For  the  Gaussian  form  |1g(P)  =  |ioexp(-p2/L2) ,  we  find(4) 


P  o 


16Vtt  L3  ( kr) '■ 


■(kL)  e 


6  -pLsm(0/2)]: 


(6) 


This  result  is  seen  to  have  the  same  wave  number  and  angular  dependence  as  that  obtained  for  the. 
square  of  the  Rayleigh  distribution  of  spherical  scatterers,  and  given  in  Eq.(13)  of  the  previous  ’ 
section,  when  we  set  L  =  2ao,  i.e.,  set  the  correlation  length  equal  to  the  rms  diameter  ot  the 
spherical  scatterers. 

To  give  some  indication  of  the  sensitivity  of  these  results  to  the  assumed  form  of  the 
correlation  function,  we  summarize  the  results  for  the  scattering  distributions  obtained  irom  the 
correlation  functions  M-e(P )  =  |ioexp(-p/L)  and  |is(p)  =  flosech  (p/L)  .  Graphs  ot  these 
expressions  are  shown  in  Fig.  3.  The  corresponding  scattering  distributions  are  found  to  be(5) 


(IpJ  Yl__Pq  1 

Vpl  2  nL3  (kr)2 

M  ^  1 

Vpl  32  L3  (kr)2 


(kL) 


[(^)2  + 1]‘ 

(^smh^sech^fe) 


(7) 


Graphs  of  the  KL  dependence  of  these  expressions,  as  well  as  that  for  the  Gaussian,  i.e.,  to(k)  - 
exp[-(KL/2)]2 ,  fE(k)=  [(KL)2  +  l]-2  and  fs(k)  =  sech2  (kKL/2)  sinh  (7tKL/2)/(7tKL/2)  ,  are 
shown  in  Fig.  4.  As  expected,  the  sharper  cutoff  for  the  Gaussian  PgIP)  is  reflected  in  a  wider 
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spectrum  for  the  transform  fc(k) .  The  negligible  difference  between  the  curves  for  fE  and  fs  is 
of  some  interest  in  its  own  right  since  it  shows  that,  at  least  in  the  present  analysis,  the  much 
maligned  cusp  at  the  origin  displayed  by  the  exponential  correlation  function  has  little  effect  on  the 
results. 

The  procedure  adopted  here,  of  assuming  a  form  for  the  correlation  function  and  then 
calculating  the  scattered  field,  is  analogous  to  that  of  assuming  a  form  for  a  refractive  index  or  a 
scattering  potential  in  deterministic  scattering.  The  reverse  procedure,  i.e.,  that  of  using  the 
scattering  data  to  determine  the  form  of  the  refractive  index  or  potential  (or  correlation  function) 
that  gives  rise  to  the  scattering,  is  referred  to  as  an  inverse  scattering  problem  and  is  a  much  more 
difficult  consideration.  For  problems  in  deterministic  scattering,  it  is  now  known  that  to  be 
successful  in  such  an  endeavor  one  must  know  much  more  about  the  scattering  data  than  is  usually 
available.  In  particular,  one  must  know  either  the  angular  distribution  of  the  scattering  at  all  angles 
at  some  one  frequency  or  know  the  scattering  at  some  one  angle  at  all  frequencies.  It  is  usually 
impracticle  to  obtain  such  information.  In  addition,  the  exact  form  of  the  refractive  index  or 
scattering  potential  may  depend  critically  on  the  scattering  data,  i.e.,  significant  changes  in  the 
predicted  form  of  the  refractive  index,  potential  or  correlation  function  may  result  from  small 
changes  in  the  scattering  data.  Finally  since  the  data  collected  may  be  expected  to  contain 
experimental  error,  it  is  usually  impossible  to  determine  any  “correct”  source  of  the  scattering. 
Fortunately,  scattering  by  random  inhomogeneities  is  usually  not  extremely  sensitive  to  the  form  ot 
the  correlation  function  and  so  a  choice  of  convenience,  frequently  the  Gaussian,  is  usually  made. 

C.  THE  INFLUENCE  OF  FLUCTUATIONS  ON  THE  DIFFRACTION  IMAGE  OF  A  FOCUSING 

SYSTEM 

The  sound  field  incident  upon  the  focusing  system  consists  of  a  coherent  part  due  to 
scattering  by  the  object  of  interest  and  an  incoherent  part  due  to  scattering  by  inhomogeneities  in 
the  medium.  When  the  amplitude  of  the  incoherent  field  becomes  comparable  to  that  of  the 
coherent  field,  the  image  can  be  assumed  to  be  no  longer  detectable.!6).  We  now  consider  the 
imaging  of  both  coherent  and  incoherent  fields. 

The  field  in  the  focal  plane  is  given  in  the  Kirchoff  approximation  by(7-8) 

i  ikr 

p{pj)  =  —  L  d'p' — ^p'Mp'’0) 

2m JS  r 


(i) 
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where  L(p')  is  the  effect  of  the  focusing  device  which  in  the  optical  case  would  be  a  lense.  The 
relation  between  r,  p’,  p  and  f  is  shown  in  Fig.  5. 

For  f  »  p,  p’  we  obtain 

r-V/'+lp-pf  (2) 

3/  +  (p=-2p.p'  +  p':)/2/ 

The  term  involving  p’2  is  usually  referred  to  as  the  Fresnel  correction.  Employing  the 
customary  procedure  of  retaining  all  terms  of  (2)  in  the  phase  in  (1)  but  setting  r  =  f  in  the 
amplitude,  which  is  less  sensitive  to  small  corrections,  we  have 

*  ikf+ikp 2  12  f 

PiPJ)  =  “  —f - 1  A*  ,2f  L(p')p(p',Q)  (3) 


If  only  the  first  term  in  the  phase  of  the  integrand  were  now  retained,  we  would  have  the 
far  field  (Fraunhoffer)  approximation  while  if  the  second  term  in  the  phase  of  the  integrand  is 
retained  we  have  the  Fresnel  approximation.  The  effect  of  the  lense  is  to  cancel  out  this  second 
term  and  yield  the  Fraunhoffer  image  in  the  focal  plane  (9).  We  thus  setrtO) 


e-ikp'2,2f 


/ 

(4) 


and  obtain 

p(p.f)  -  '!/,l  (5) 

The  field  in  the  focal  plane,  p(p,f)  ,  has  now  been  expressed  as  an  integral  over  the  field  in  the 
plane  of  the  receiver  p(p,0) . 

The  field  in  the  receiver  plane,  p(p,0),  may  be  considered  to  be  composed  of  a  coherent 
part  that  represents  the  information  of  interest  and  a  random  part  that  is  due  to  scattering  by 
inhomogeneities  in  the  medium.  As  in  Section  B,  we  treat  the  randomness  of  the  pressure  wave  by 
calculating  I p ( p . f ) 1 2  ,  a  quantity  proportional  to  the  intensity  of  the  wave,  and  then  introduce  a 
correlation  function  to  describe  the  fluctuations  that  are  now  associated  with  the  sound  field  at  tire 
receiver  rather  than  with  fluctuations  in  the  properties  of  the  medium. 

From  (5)  we  have 
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(|/>(p,/)|2)  = 


2  mf 


Is  M 


-*p«(p'-p  ")// 


^(|p'-p"l) 


where,  following  the  procedure  used  in  part  B,  we  have  set 


^(Ip,-p1)=(p*(p'-0)p(p",o)) 


(6) 


(7) 


As  noted  in  Section  B,  this  choice  for  a  correlation  function  implies  an  assumption  that  the 
fluctuating  quantity,  here  the  sound  field,  is  both  homogeneous  and  isotropic. 

The  coherent  field.  The  above  result  is  immediately  specialized  to  the  case  of  a  coherent 
field  by  setting  |i(|p'-p"|)  =  p02  where  po2  is  proportional  to  the  intensity  of  the  incident  wave. 
In  this  limit  there  is  no  coupling  between  p'  and  p"  and  we  have 


(Mp>/)1  )_  i .[.  d2  ,w// 

(p0kA/2nf)2  A2  U* 


A  =  Area  of  Receiver 


ijkpalj)- 

kpa / 2/ 

sin(kxh/2f) 
~kxh  /  2/ 


,  circular  receiver,  A  =  rra2 


sm(kyh!2f) 

icyhJYf 


square  receiver,  A  =  h2 


✓ 


(8) 


where  Ji  refers  to  the  Bessel  function  of  first  order.  Both  of  these  expressions  have  a  form  that 
peaks  at  unity  at  the  center  of  the  receiver.  The  results  for  the  two  geometries  are  compared  in  Fig. 
6  in  which  (|pC0hl2)/pc2>  with  pc  =  p0kA/27tf,  is  plotted  as  a  function  of  distance  from  the  center 
for  both  a  circular  and  a  square  receiver.  The  square  case  is  considered  because,  as  summarized 
below,  it  is  possible  to  obtain  a  simple  complete  description  of  the  incoherent  field  for  this 
geometry(ll>. 

For  the  circular  receiver,  we  first  set  k  =  2itl\  and  write  the  argument  of  the  Bessel 
function  as 

kpa  .  a  a  p 

/  a 


(9) 
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The  choices  A.  =  0.3  meter  (which  follows  from  a  sound  velocity  of  300  meters/sec.  and  a 
frequency  of  1000  Hz.)  as  well  as  a  =  0.5  meter  and  f  =  1  meter,  have  been  used.  The  result  is 
given  by  the  curve  marked  “1”  in  Fig.  6. 

For  the  square  receiver  we  consider  distance  from  the  center  along  a  diagonal  x  =  y 
indicated  by  curve  “2”  in  Fig.  6,  as  well  as  along  a  coordinate  line  (y  =  0,  0  <  x  <  h/2).  The 
latter  curve  is  indicated  by  “3”  in  Fig.  6.  (Since  on  the  diagonal  a  radial  distance  R  from  the 
center  corresponds  to  a  value  of  x  equal  to  R  cos  k/4  =  R/V2,  for  the  graph  marked  “2”  x  has 
been  replaced  by  x/V2  to  account  for  this  scaling.)  As  expected,  near  the  center  curves  2  and  3 
yield  equal  values  for  the  coherent  field  intensity.  The  correspondence  between  circular  and  square 
receivers  is  also  found  to  be  quite  close  as  should  also  be  expected.  For  numerical  values  in  the 
rectangular  case  we  have  used  h/2  =  0.5  meter  with  the  other  parameters  remaining  the  same  as  for 
the  circular  case. 

Incoherent  field.  Here  we  must  make  a  choice  for  the  correlation  function.  As  noted 
above,  the  choice  is  a  balance  between  physical  realism  and  analytical  tractability  and  the  Gaussian 
is  frequently  chosen.  We  adopt  this  form  here  and  set 


MIp'-pI  )  =  pU 


-|p'-p"|:  /L: 


(10) 

/ 


where  pp2  is  proportional  to  the  mean  square  intensity  of  the  fluctuating  field  and  L  is  the 
correlation  length,  i.e.  the  distance  over  which  the  fluctuations  in  the  field  display  some  coherence. 
We  are  thus  introducing  a  new  length  into  the  considerations  and  can  examine  results  as  a  function 
of  the  ratio  of  the  receiver  size,  a  or  h,  to  the  correlation  length  L.  For  arbitrary  ratios  of 
receiver  size  to  correlation  length  it  is  found  that  the  results  can  be  obtained  in  a  simple,  readily 
interpretable  form  only  for  the  rectangular  receiver.  For  the  circular  case  it  appears  that  only  the 
large  receiver  limit  a/L— ><»  can  be  treated  in  this  manner^).  The  results  are  of  considerable 
interest  since  they  show  that  the  fluctuations  in  the  image  are  diminished  by  increasing  the  receiver 
size  until  it  is  about  twice  the  correlation  length  L.  This  result  is  analogous  to  the  optical  effect  that 
the  twinkling  of  stars  appears  greater  to  the  naked  eye  than  to  a  larger  optical  receiver. 

For  the  square  case  we  write 

p  •  (p'  -  p")  =  x{x'  -  x")  +  y(/  -y")  UU 

and  express  (6)  in  the  form 


' kh ^ 


KoMff)  =  P?[^\  M{h/L7x)M(hl 


(12) 
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where 


M(h/L,r,)  = 


1  rhn 


cjx"e-*x(x‘-x‘)i  f-v-x")1  nl 


(13) 


and  2yx  =  kxL/f.  The  factor  M(h/L,yy)  has  the  same  form  as  M(h/L,yx)  except  that  x  is 
replaced  by  y  .  These  quantities  have  been  defined  so  that  they  reduce  to  unity  at  the  center  of  a 
very  large  receiver,  i.e.,  M(°°,0)  =  1.  (It  should  be  noted  that  results  for  a  rectangular  receiver 
may  be  obtained  by  merely  modifying  the  limits  on  the  integral  for  either  M(h/L,yx)  or 
M(h/L,yy). 

The  details  of  the  integration  are  summarized  in  Appendix  A.  One  obtains 


M{h/L,yx)  =  2{if{±P-Q) 


(14) 


where 

P  =  \4ne~rx  Re{F(h/  L,yx)} 

Q  =  \[  1  -  e~(h'L) 2  cos(2 yxh!  L)]  +  \4%e~rl  lm{F(h/ L,yx)}  (15) 

in  which 

'  F(h  /  L,yx)  =  erf(j  +  iyx )  -  eif{iyx )  (16) 

and 


e,f{z)  =  r  I'  dte-2 

J  V  '  SK  Jo 

(17) 

Results  of  the  same  form  occur  for  the  y  dependence  with  yy  =  kyL/2f  being  used  in  place  of 

Yx- 

To  obtain  some  understanding  of  this  result,  we  now  examine  various  special  cases.  It  is 
convenient  to  use  the  coefficient  in  (12)  as  a  reference  value  and  set 


12 


so  that 


(\p,ncoh(pjf) 

Pr 


M(h/L,Yx)M(h/L,yy) 


(18) 


(19) 


1)  Field  at  the  center  of  the  receiver  (x  —  y  =  0)  as  a  function  of  h/L  In  this  case,  yx 
=  Yy  =  0  and  from  (19)  we  have 


(k„.(o./)f) 

Pr 


[M(h/L,  0)]2 


(20) 


The  results  in  (14)  to  (16)  yield 


M{h/L,  0)  =  Vn^j 


+ 


e-vuLf 


,  (21) 


as  obtained  by  Chernovd3).  A  graph  of  (20)  is  shown  as  the  solid  curve  in  Fig.  7.  A  graph  of 
exp[-(h/L)2]  has  been  included  for  comparison.  The  fluctuations  are  seen  to  attain  their  maximum 
value  in  the  limit  of  small  receiver  size  to  correlation  length  (h/L  «  1)  and  to  decrease  rapidly  as 
h/L  increases  until  a  ratio  of  h/L  equal  to  about  two  is  obtained,  after  which  they  decrease  more 
slowly.  As  noted  above,  this  decrease  in  the  effect  of  fluctuations  is  well  known  in  astronomical 
observations  where  the  scintillation  of  stars  is  found  to  be  much  larger  to  the  naked  eye  than  when 
observed  with  large  optical  receivers. 

2)  Noncentral  points  To  examine  the  corresponding  situation  at  various  non  central 
receiver  locations,  we  write  pincoh(p1f)  =  Pincoh(x,y,f)  and  consider,  in  Fig.  8,  the  value  of 
<  IPincoh(x,x,f)|2  >/pR2,  i.e.  the  field  at  locations  along  the  diagonal  x  =  y,  for  x  =  O.lh,  0.3h 
and  0.5h  . 

In  obtaining  these  results  we  have  had  to  make  a  choice  for  the  parameters  involved  in  the 
definition  of  2yx  =  kxL/f.  Setting  k  =  2ro9i,  we  first  write 


_  _  L  L  x  h 

2r-=2*I7*I 


(22) 
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The  choices  X  -  .3  meters  and  f  =  0.2  meters  have  again  been  used  as  well  as  L=  0.1  meter. 
No  attempt  at  a  thorough  parameter  study  has  been  made  at  this  stage.  Although  the  result  in  Fig. 
8  predicts  a  more  rapid  decrease  in  the  fluctuations  as  one  moves  away  from  the  center,  there  is 
probably  no  advantage  to  be  gained  from  this  decrease  since  the  coherent  field  itself  falls  off  as 
one  moves  away  from  the  center  (cf.  Fig.  6). 

3)  Spatial  Dependence  for  fixed  h/L  We  first  consider  the  case  of  large  receiver  to 
correlation  length  ratio  (h/L  »  1).  This  limit  is  easily  treated  by  merely  setting  h/L  =  oo  in  the 
upper  limit  of  the  integrals  that  define  P  and  Q  as  given  in  Appendix  A.  We  then  obtain 


M(h/ L,y) 


hi L»y  (=/,,yv) 


(23) 


and  find 

(I  PU^yjf) 

p\  h2 

This  result  has  the  same  functional  form  as  that  given  by  IshimaruO4)  for  the  large  circular  receiver 
with  x2  +  y2  =  r2.  As  noted  by  ChernovO2',  far  from  the  edge  of  the  receiver,  where  shape 
effects  should  be  negligible,  results  for  both  square  and  circular  receivers  can  be  expected  to  be 
quite  similar. 

For  h/L  =1,2  and  4  and  field  points  on  the  diagonal  (i.e.,  y  =  x),  we  obtain  the  results 
shown  in  Fig.  9.  The  decrease  in  the  magnitude  of  pincoh  associated  with  increasing  values  of 
h/L  is  consistent  with  the  earlier  results  for  dependence  of  fluctuations  on  receiver  size  shown  in 
Fig.  7. 
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Appendix  A  Evaluation  of  M(h/L,yx) 


From  (12)  of  Section  C  we  have 


dx'l 


tf  -ikx(x'-x")/ f-(x'-x")2  /  l} 


<fxV* 


(A.  1) 


Consider  this  integral  as  being  taken  over  a  square  of  area  h2  in  an  x'-x"  plane.  Separating 
variables  in  the  integrand  by  introducing  u  =  (x'  -  x")  and  v  =  (x'  +  x")/2,  a  transformation  for 
which  dx'dx"  =  dudv,  we  integrate  over  the  same  area  in  the  u-v  plane  by  evaluating 


h2M(hl  L,yx)=f  du\Vl+U)ri  dve-i{hc,f)u-{u,L)1  +  f  duf'~U>12  dye" 
V  '  x>  J-/,  J-(/j+u)/2  JO  J-(h-u)l2 


-i(kxf  f)u-(uJ  L)~ 


(A.2) 


The  area  is  shown  in  Fig.  10  where  two  examples  of  the  area  strips  that  are  summed  in  this 
integration  process  are  indicated.  When  the  v  integration  is  carried  out  we  find 


/ 

(A.3) 


where 

P  —  |  dE,e~^  cos2^;  =  Rejjo  d^e"1’ e2,y<; | 

{2  =  Jo  *d<i;<r^£cos2y£  =  Re|jo  Ldge~r  £e2,}*J  (A-4) 

with  ^  =  u/L  and  2y  =  kxL/f.  These  integrals  may  be  expressed  in  terms  of  the  complex  error 
function.  The  results  are  given  in  (14)  -  (17). 

For  h/L  »  1,  we  replace  the  upper  limit  on  each  integral  in  (A.4)  by  °°  and  have 


F  =  Jo  di;e  ^  cos2)^  =  j-sfne  r 

Q  =  j0  d%e~*  l; cos 2y£  =  Jo  dfy~V  %  1-^#-  +  ...  «^--y2+... 


Since  P  is  multiplied  by  h/L  in  (A.3),  the  contribution  from  Q  is  negligible  unless  y»l. 


15 


References 

1 .  L.  Amundsen,  B.  Amtsen  and  R.  Mittet,  “Depth  Imaging  of  Offset  Vertical  Seismic  Profile 
Data,”  Geophysical  Prospecting,  41,  1009-1031,  (1993). 

2.  If  spatial  variations  in  the  density,  rather  than  the  bulk  modulus,  were  considered,  we 
would  obtain  the  more  complicated  wave  equation 

d2P 

-r ~  ~  ( B/p)V2P  +  (B/ p)V(ln  p)  •  VP  =  0 

at 

At  this  preliminary  stage  in  the  analysis  where  only  a  qualitative  understanding  of  the 
problem  is  being  developed,  it  seems  appropriate  that  this  more  complicated  wave  equation, 
along  with  the  shear  wave  considerations,  should  be  neglected. 

3.  A.  Erdelyi,  ed.  Tables  of  Integral  Transforms,  McGraw  Hill,  New  York,  (1954),  Vol.  2, 
p.  29,  No.  10. 

4.  Ref.  3,  Vol.  1,  p.  73,  No.  19. 

5.  Ref.  3,  Vol.  1,  p.  152,  No.  14  and  p.  89,  No.  12. 

6.  A.  Ishimaru,  Wave  Propagation  in  Random  Media,  Academic  Press,  New  York  (1978), 
Vol.  2,  p.  305. 

7.  Ref.  6,  p.  302. 

8.  M.  Born  and  E.  Wolf,  Principles  of  Optics,  3rd.  (Revised)  Ed.,  Pergamon  Press  (1965), 
New  York,  pp.  382  -  385. 

9.  Ref. 8,  pp.  384  -  385. 


10. 


Ref.  6,  p.  302. 


16 


11.  L.  A.  Chernov,  Wave  Propagation  in  a  Random  Medium,  Dover  (1960),  New  York,  p. 
137. 


12.  Ref.  6,  p.304. 

13.  Ref.l,  p.  153. 

14.  Ref.  6,  p.  304. 


/ 


P'i-6 


F/j./a 


